Shrinking Magnetic Vortices in V 3 Si Due to Delocalized Quasiparticle Core States: 
Confirmation of the Microscopic Theory for Interacting Vortices 
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We report muon spin rotation measurements on the conventional type-II superconductor VsSi that 
provide clear evidence for changes to the inner structure of a vortex due to the derealization of bound 
quasiparticle core states. The experimental findings described here confirm a key prediction of recent 
microscopic theories describing interacting vortices. The effects of vortex- vortex interactions on the 
magnetic and electronic structure of the vortex state are of crucial importance to the interpretation 
of experiments on both conventional and exotic superconductors in an applied magnetic field. 

PACS numbers: 74.20.Fg, 74.25. Qt, 74.70.Ad, 76.75.+i 



In 1964, a breakthrough paper by Caroli, de Gennes 
and Matricon p| showed that in the framework of the mi- 
croscopic theory, quasiparticles (QPs) bound to an iso- 
lated vortex of a conventional s-wave type-II supercon- 
ductor occupy discrete energy levels. Twenty-five years 
later, localized vortex core states were observed for the 
first time in NbSe2 by scanning tunneling microscopy 
(STM) @. Our understanding of the vortex state in 
type-II superconductors has accordingly progressed from 
Abrikosov's initial prediction Q based on the macro- 
scopic Ginzburg-Landau (GL) theory Q, to current the- 
ories describing the electronic structure of magnetic vor- 
tices on a microscopic level. However, it is only in re- 
cent years that predictions have emerged from the mi- 
croscopic theory on the effects of vortex-vortex interac- 
tions. In analogy with bringing atoms close together to 
form a conducting solid, increasing the vortex density by 
applying a stronger magnetic field H enhances the over- 
lap of bound state wave functions of neighboring vor- 
tices, resulting in the formation of energy bands that 
allow the intervortex transfer of QPs H, H S 13 • This 
is expected to strongly influence experiments on conven- 
tional superconductors that are sensitive to quasiparticle 
excitations, such as specific heat and thermal transport, 
and to have profound effects on the magnetic structure 
of the vortex state. However, understanding the poten- 
tial interplay between vortices and quasiparticles is also 
of crucial importance in the study of high-temperature 
superconductors. In these and other exotic superconduc- 
tors, comparatively little is known about the structure 
of the vortex state and its effect on experiments in large 
magnetic fields. It is therefore essential to have a solid 
understanding of the behavior of interacting vortices in 
conventional superconductors, and to establish the con- 
nections with quasiparticle properties. 

The effect of delocalized QP core states on the spa- 
tial variation of the pair potential A(r) at a vortex site 



has been considered in the framework of the quasiclas- 
sical Eilenberger theory 0, H, @- These calculations 
show that the effect of the intervortex transfer of QPs on 
A(r) leads to a reduction of the size of the vortex cores 
with increasing H. Such shrinking of the vortex cores 
has in fact been observed by muon spin rotation (/iSR) 
H3 [H [H [H Q E and STM HI , and proposed 
as a possible explanation for the anomalous low-field be- 
havior observed of the specific heat in conventional su- 
perconductors 0, 0. However, experimentally there 
has been inadequate evidence in support of a causal re- 
lationship between the size of the vortex cores and the 
derealization of bound QP core states. 

An important finding has come by way of recent low- 
temperature thermal conductivity measurements, which 
are sensitive only to extended or delocalized electronic 
excitations. These studies have revealed the existence of 
highly delocalized QPs down to low magnetic fields in the 
vortex state of LuNi 2 B 2 C [U, YBa 2 Cu 3 6 .9 HI, NbSe 2 
[22j| and MgB 2 :23| . In the extreme gap anisotropy super- 
conductors LuNi 2 B 2 C and YBa 2 Cu306.g, the dominant 
contribution to the thermal conductivity is believed to 
be a field-induced Doppler shift of the QP spectrum out- 
side the vortex cores [24j . Most surprising are the results 
for NbSe 2 , long believed to be representative of a simple 
conventional s-wave superconductor. The high degree of 
QP derealization in NbSe 2 appears to arise from two- 
gap superconductivity [13, as is the case for MgB 2 [23| . 
In these superconductors, the heat at low fields is car- 
ried by the QP excitations associated with the smaller 
of the two energy gaps. This smaller gap is a possible 
source of the very large low-field value of the core size 
observed in NbSe 2 by /xSR [13 and in MgB 2 by STM 
H3- On the other hand, the large core size observed at 
low fields in borocarbide 0, 0] and cuprate 0, 0, Il3| 
superconductors is not completely understood. 

The conventional superconductors VaSi and Nb pro- 



vide a unique opportunity to experimentally observe the 
effects of intervortcx QP transfer on the inner structure 
of a vortex. In contrast to the abovementioned super- 
conductors, thermal conductivity measurements p3. l26| 
reveal no appreciable derealization of QP states at low 
H — as expected when the vortices are nearly isolated. 
VjSi is particularly suitable, since it has a simple cubic 
crystal structure, a relatively high superconducting tran- 
sition temperature (T c = 17 K) and a large upper criti- 
cal field (H c2 ~ 185 kOe). To determine the size of the 
vortex cores we measured the internal magnetic field dis- 
tribution in VaSi by /iSR at the Tri-University Meson 
Facility (TFJUMF), Vancouver, Canada. The experi- 
ment was performed by implanting spin-polarized muons, 
which stop randomly in the sample on the length scale 
(~ 10 2 -10 3 A) of the vortex lattice (VL). The magnetic 
moment of the muon precesses about the local magnetic 
field B with a Larmor frequency given by v = j^B, where 
7 M = 0.0852 /is _1 G _1 is the gyromagnetic ratio of the 
muon. The spatial distribution of B is determined by 
measuring the time evolution of the muon spin polariza- 
tion via the anisotropic distribution of decay positrons 

M 

Small-angle neutron scattering |2g and recent STM 
[29| images of VjSi show that for a field applied along 
the fourfold [001] axis, the VL undergoes a transition 
from hexagonal to square symmetry. Kogan et al. |30| 
have developed a phenomenological London model that 
attributes this transition to nonlocality of the relation 
between the supercurrent density j(r) and the vector po- 
tential A in a region around the vortex core. This estab- 
lishes a connection between the vortex structure and the 
anisotropy of the Fermi surface, which in VaSi produces 
fourfold symmetry near the vortex cores. We note that 
a more general model has recently been developed for 
fourfold symmetric superconductors, which incorporates 
the effects of both Fermi surface and gap anisotropy on 
the vortex structure At fields below 7.5 kOe where 
the intervortex spacing is large, the isotropic magnetic 
repulsion of vortices yields a hexagonal lattice. With in- 
creasing field there is an increased overlap of the fourfold 
symmetric regions of the individual vortices, and the VL 
evolves into a square at 40 kOe. 

The Kogan model uses an arbitrary cutoff factor to 
overcome the logarithmic divergence of the internal mag- 
netic field B(r) at the vortex sites. We find that the /zSR 
data are fit much better using a GL analog of the Kogan 
model, which properly accounts for the finite size of the 
vortex cores [32J . The spatial profile of the periodic local 
magnetic field in the GL formalism is given by 



B(r) = B Q (l-b 4 ) 
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u K\{u) 



X 2 G 2 + \ 4 (n xxyy G 4 + dG 2 G 2 ) 



(1) 

Here b — B/B c2 is the reduced field, B is the average 
internal field, G are the reciprocal lattice vectors, K\{u) 
is a modified Bessel function, u 2 = 2£ 2 G 2 (1 + b 2 )[l — 
26(1 — b) 2 ], £ is the GL coherence length, and n xxyy and 
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FIG. 1: Fourier transform (FT) of the measured muon spin 
precession signal (green circles) in V3S1 at T = 3.8 K and 
H = 50 kOe. The FT provides an approximate illustration 
of the internal magnetic field distribution, but is broadened 
by the apodization procedure used to smooth out the ringing 
effects of the finite time window (0-3 fjs) and the noise from 
reduced counts at later times. Additional broadening due to 
VL disorder and nuclear magnetic dipoles is accounted for by 
a Gaussian broadening width of <r~ 1.0 fJS ■ The solid red 
curve is the FT of the fit in the time domain, assuming the 
field profile B(r) given by Eq. The inset is a contour plot 
of the function B(r) obtained from the fit. 



d are dimensionless parameters arising from the nonlo- 
cal corrections. The quartic n xxyy term is an isotropic 
correction, whereas the biquadratic d term controls the 
fourfold anisotropy. 

Figure 1 shows Fourier transforms of both the mea- 
sured muon-spin precession signal and the fit to the data 
in the time domain at 50 kOe applied parallel to [001]. 
In agreement with the imaging experiments, we obtained 
excellent fits to a square and hexagonal VL above 40 kOe 
and below 7.5 kOe, respectively, and to a rhombic unit 
cell with an apex angle (3 w 60° — 90° at intermediate 
fields. These fits and those to Kogan's model also yielded 
the following results: (i) For all H, n xxyy « 0. This is con- 
sistent with an analogous model developed by Affleck, 
Franz and Amin |33j ] , who found that the quartic correc- 
tion makes a negligible contribution, (ii) Below 7.5 kOe, 
excellent fits were obtained with d = 0. This indicates 
that the VL is not strongly tied to the underlying crystal 
lattice, as was determined by STM. (iii) Above 7.5 kOe, 
the orientation of the vortex cores determined by our 
analysis (see Fig. 1 inset) is that which is expected for 
close packing of square vortices. 

The field dependence of A, £, d and (3 are shown in 
Fig. 2. At low fields where the lattice is hexagonal, 
A « 1060 A and £ « 42 A, which are consistent with pre- 
viously determined values of these parameters. Above 
7.5 kOe there is a slight nonphysical increase in the fit- 
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FIG. 2: The magnetic field dependence of the fit parame- 
ters at T = 3.8 K. (a) The magnetic penetration depth A. (b) 
The coherence length f. (c) The apex angle f3 (purple cir- 
cles) and the anisotropy parameter d (green triangles). The 
dashed vertical lines indicate the field range over which the 
VL continuously evolves from hexagonal to square symme- 
try. At fields immediately above 7.5 kOe, d and (5 are poorly 
determined, because the VL is only slightly distorted from 
hexagonal symmetry. 



ted values of A and £, reflecting the beginning of the 
gradual transition to a square VL. At higher fields £ de- 
creases, and continues to do so when the VL completes 
its transition to a square at 40 kOe. 

Gygi and Schluter [34| have shown that the pair poten- 
tial A(r), which in general is a numerical function, varies 
on two length scales, £i and £2- The first is defined from 
the slope of A(r) near the center of the vortex core 

AM 

6 = A / lim -LI (2) 

r^O r 

(Ao is the BCS superconducting energy gap), whereas £2 
is the length scale over which A(r) rises to its asymptotic 



value Ao- It is important to note that while at low tem- 
peratures £2 is close to the value of the coherence length 
derived from H c2 (i.e. ~ 42 A), £1 is considerably shorter 
[see Fig. 3(b)]. Thus, we can understand the reduction 
of £ above 7.5 kOe as a rapid variation of A(r) at the 
vortex center, and not an increase of H c i- Theoretically, 
the rapid increase of A(r) near r — is accompanied 
by a similar rise of the supercurrent density j(r) in the 
same region 0, H, 13 ■ Consequently, we can define an 
effective vortex-core size ro as the distance from the core 
centre (r = 0) to the location where j(r) reaches its max- 
imum value, measured along the line connecting nearest- 
neighbor vortices. Solutions of the microscopic theory 
show that £1 and ro exhibit the same qualitative behavior 
as functions of H and T. From previous /iSR studies of 
the Kramer-Pesch effect, it is known that the fitted value 
of £ in Eq. and ro exhibit similar behavior, and hence 
£ reflects a particular sensitivity to £1. The parameter ro 
can be obatined in a nearly model-independent way us- 
ing the Maxwell relation, j(r) = |V x B(r)|, where B{r) is 
obtained from fitting the /iSR time spectrum. Because ro 
is not a fit parameter, the details of the theoretical model 
for B(r) are not very important in this procedure. What 
is essential is that the model for B(r) yields excellent fits, 
as was the case using Eq. (QJ. 

In Fig. 3(a) we compare the field dependence of ro in 
VsSi to the electronic thermal conductivity n e measured 
previously j^. For the first time we see the expected 
field dependence of the core size in a single-gap conven- 
tional s-wave superconductor. In contrast to NbSe2 01 
and MgB2 ;25j , the low-field value of ro is consistent with 
the coherence length calculated from H C 2- At low fields, 
where only a modest overlap of the QP core states of 
neighboring vortices is expected, both ro and K e exhibit 
a weak dependence on H. The change in ro over this 
range of H is primarily due to the superposition of the 
j(r) profiles of nearest-neighbor vortices. This is indi- 
cated by the green circles in Fig. 3. Above 7.5 kOe the 
core size shrinks more rapidly than that expected from 
the superposition of j(r) profiles, and is accompanied by 
a simultaneous increase in the electronic thermal conduc- 
tivity. Together these observations signify a change in the 
slope of A(r) at the vortex center, due to an increased 
overlap of the bound state wave functions of adjacent 
vortices. As shown in Fig. 3(a), the transformation to a 
square vortex lattice begins at the field where the der- 
ealization of QPs becomes significant. In other words, it 
is the increased strength of the vortex- vortex interactions 
that drives the symmetry change of the lattice to reflect 
the fourfold symmetry of the individual vortices. 

The experimental results presented here for VaSi con- 
firm one of the key predictions of theoretical works that 
advocate the importance of intervortex QP transfer to 
the detailed structure of the VL. This important detail 
should be a consideration in the interpretation of any 
experiment on a type-II superconductor in an applied 
magnetic field. 
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FIG. 3: (a) The magnetic field dependence of the vortex core 
size ro measured by /iSR (blue circles) at T — 3.8 K, and 
the electronic thermal conductivity K e /T (red squares) from 
Ref. [2^| extrapolated to T — > (and normalized to the value 
kn/T at H C 2). The green circles indicate the reduction of 
ro due to the superposition of j(r) profiles from individual 
vortices, calculated assuming the low-field value of £. The 
data are plotted as a function of 1/H 1 ^ 2 , which is proportional 
to the intervortex spacing L = (&o/H sin 0) 1/72 , where <E>o is 
the magnetic flux quantum. The dashed vertical lines indicate 
the field range over which the VL undergoes a continuous 
hexagonal-square transition. Schematic of the pair potential 
A(r) as a function of distance from the vortex core centre at 
high (b) and low (c) magnetic field. 
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